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We present a numerical model for a two dimensional (2D) 
granular assembly, falling in a rectangular container when the 
bottom is removed. We observe the occurrence of cracks split- 
ting the initial pile into pieces, like in experiments. We study 
in detail various mechanisms connected to the 'discontinuous 
decompaction' of this granular material. In particular, we 
focus on the history of one single long range crack, from its 
origin at one side wall, until it breaks the assembly into two 
pieces. This event is correlated to an increase in the number 
of collisions, i.e. strong pressure, and to a momentum wave 
originated by one particle. Eventually, strong friction reduces 
the falling velocity such that the crack may open below the 
slow, high pressure 'dynamic arch'. Furthermore, we report 
the presence of large, organized structures of the particles' 
angular velocities in the dense parts of the granulate when 
the number of collisions is large. 
PACS: 46.10.+Z, 05.60+w, 47.20.-k 



I. INTRODUCTION 

The flow behavior of granular media in hoppers, pipes or 
chutes has received increasing interest during the last years. 
For a review concerning the physics of granular materials, see 
and refs. therein. 
The complete dynamical description of gravity driven flows 
is an open problem and for geometries like hoppers or verti- 
cal pipes, several basic phenomena, are still unexplained. In 
hoppers intermittent clogging due to vault effects j|, density 
waves in the bulk @ or 1// noise of the outlet pressure, have 
been reported from experiments || . In a vertical pipe geom- 
etry, numerical simulations on model systems with periodic 
boundary conditions and also analytical studies Q show 
density waves. But so far, experimental evidence of density 
waves is only found in situations where pneumatic effects, i.e. 
gas-particle interactions, are important 

In the rapid flow regime, kinetic theories jnj,[nj describe 
the behavior of the system introducing the granular temper- 
ature as a measure of the velocity fluctuations. In quasistatic 
situations, also arching effects and particle geometry get to 
be important. Under those conditions, kinetic theories and 
also continuum soil mechanics do not completely describe 
all the phenomena observed. The complicated particle-wall 
and particle-particle interactions |l2"|-[l4"[ , possible formation 
of stress chains in the granulate and also stress fluctuations 



l. r )||l(j| | undoubtedly require further experimental, theoretical 
and simulation work. 

Following recent research on model granular systems 
[p^jil]] we focus here on the problem of a 2D pile made up of 
rather large spheres enclosed in a rectangular container with 
transparent front and back walls, separated by slightly more 
than one particle diameter. Recent observations of approxi- 
mately V-shaped microcracks in vertically vibrated sand-piles 
fll7| were complemented by recent experiments and simula- 
tions of the discontinuous decompaction of a falling sandpile 
p2~[ . From experiments Duran et al. ]z| find the following ba- 
sic features: In a system with polished laterals walls, cracks 
are unlikely to appear during the fall, i.e. the pile will acceler- 
ate continuously, with an acceleration value depending on the 
aspect ratio of the pile and on the friction with the walls. In a 
system with rather poorly polished walls (a surface roughness 
of more than 1 /jm in size) , cracks occur frequently. A crack in 
the lower part of the pile grows, whereas a crack in the upper 
part is unstable and will eventually close. These experimen- 
tal findings can be understood from a continuum approach 
based on a dynamical adaptation of Janssen's model [^,p2[. 
Nevertheless, not all of the discontinuous phenomena, such as 
the reasons for the cracks, can be explained by such a contin- 
uum model. In previous works, numerical simulations were 
used to parallel the experiments and to analyze the falling 
pile for different material's parameters pipi) ]. Though the 
simulations are dynamic, in contrast to an experiment which 
starts from a static situation, a reasonable phenomenological 
agreement was found. Furthermore, simulations were able to 
correlate the existence of long range cracks to strong local 
pressure on the walls. The increase in pressure was found to 



In this work we use 
5| and investigate in 



be about one order of magnitude 
the numerical model of Refs. J24|j22 
detail how a crack occurs. In particular, we follow one spe- 
cific crack and try to extract the generic features that could 
be relevant for a more involved theoretical description of the 
behavior of granular materials. 

We briefly discuss the simulation method in section [n| and 
present our results in section III which consists of three main 
parts: Firstly, we present stick-slip behav ior and a local or- 
ganization of the spins in subsection III A . Secondly, we de- 



scribe in how far the number of c ollisio ns, the kinetic energy 
and the pressure are connected in [II B and finally we present 
long range organization and momentum waves in subsection 
1IIC. We summarize and conclude in section £v|. 
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II. THE SIMULATION METHOD 

Our simulation model is an event driven (ED) method 
[ p4||2^ , p5| -fjj| based upon the following considerations: Parti- 
cles undergo a parabolic flight in the gravitational field until 
an event occurs. An event may be the collision of two particles 
or the collision of one particle with a wall. Particles are hard 
spheres and interact instantaneously; dissipation and friction 
are only active on contact. Thus we calculate the momentum 
change using a model that is consistent with experimental 
measurements [ p8| . From the change of momentum we com- 
pute the particles' velocities after a contact from the velocities 
just before contact. We account for energy loss in normal di- 
rection, as for example permanent deformations, introducing 
the coefficient of normal restitution, e. The roughness of sur- 
faces and the connected energy dissipation, is described by 
the coefficient of friction, fj,, and the coefficient of maximum 
tangential restitution /3q. For interactions between particles 
and walls we use an index w, e.g. Due to the instanta- 

neous contacts, i.e. the zero contact time, one may observe 
the so-called 'inelastic collapse' in the case of strong dissipa- 
tion. For a discussion of this effect see McNamara and Young 
[p9| and refs. therein. Despite this problem, we use the hard 
sphere ED simulations as one possible approach, since also 
the widely used soft sphere molecular dynamics (MD) simu- 
lations may lead to complications like the 'detachment-effect' 
Jjo|-|32[ or the so-called 'brake-failure' for rapid flow along 
rough walls |Q . Recent simulations of the model system, de- 
scribed in this study, using an alternative simulation method, 
the so called 'contact dynamics' (CD) , also lead to cracks 
[ pij . The CD method has a fixed time-step, in contrast to 
ED where the time-step is determined by the time of the next 
event. 

From the momentum conservation laws in linear and angu- 
lar direction, from energy conservation, and from Coulomb's 
law we get the change of linear momentum of particle 1 as a 
function of e, 11, and Po Pq]: 



AP = - mi2 (l + e )t£ 
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with the reduced mass m\2 = mim2/(mi +7712). For particle- 
wall interaction, we set rri2 = 00 such that ui\2 = m 1 . (n) 
and (t) indicate the normal and tangential components of the 
relative velocity of the contact points 
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with Vi and being the linear and angular velocities of par- 
ticle i just before collision, di is the diameter of particle i 
and the unit vector in normal direction is here definded as 
n = (fx — f*2)/|fi — F2I. Paralleling e, the (constant) coefficient 
of normal restitution we have (3, the coefficient of tangential 
restitution 



with the collision angle n/2 < 7 < it |25|. Here, we simplifed 
the tangential contacts in the sense that exclusively Coulomb- 
type interactions, i.e. AP'*' is limited by /iAP' n \ or sticking 
contacts with the maximum tangential resitution /3o are al- 
lowed. Sticking corresponds thus to the case of low tangential 
velocity whilst the Coulomb case corresponds to sliding, i.e. 
a comparatively large tangential velocity. For a d etailed dis- 
cussion of the interaction model used see Refs. $24 



III. RESULTS 

Since we are interested in the falling motion of a compact 
array of particles, we first prepare a convenient initial condi- 
tion. Here, we use N — 1562 particles of diameter d = 1 mm 
in a box of width L — 20. 2d, and let them relax for a time t r 
under elastic and smooth conditions until the density and en- 
ergy profiles do not change any longer. The choice of L is quite 
arbitrary, however, we wanted to start with a triangular lat- 
tice with a lattice constant of d(l + A) and A to be small but 
larger than zero, i.e. A = 0.01. We tested several height to 
width ratios S = H/L of the system and found e.g. the same 
behavior as in experiments, i.e. the larger S, the stronger are 
the effects discussed in the following. The average velocity of 
the resulting initial condition is v = V< v 2 > « 0.05m/s. 




j3 = min [/3o,/?i 



(3) 



/3o is the coefficient of maximum tangential restitution, — 1 < 
/3o < 1, and accounts for energy conservation and for the elas- 
ticity of the material [Q. /3i is determined using Coulomb's 
law such that for solid spheres /3i = —1 — (7/2)/i(l + e) cot 7 
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t = 0.048 s 0.056 s 0.064 s 



FIG. 1. (a) Snapshots of a typical experiment (left) and a 
simulation (right) with N =1562 particles, at time t — 0.06 s, 
in a pipe of width L/d = 20.2. We used here e = 0.96, e w = 
0.92, fi = 0.5, (Mb = 1.0, and (3 = (3 0w = 0.2. (b) Snapshots 
from a simulation with almost the same parameters as in (a), 
but here e = 0.90, e w — 0.90, and (i w = 0.5. The greyscale 
indicates the number of collisions per particle per millisec- 
ond. Black and white correspond to no collision or more than 
ten collisions respectively, (c) Here only the selected parti- 
cles #369, #370, #371, #390, #409, and #428 from (b) are 
plotted to illustrate their motion. The vertical line indicates 
the right wall. 



The kinetic energy connected to v is comparable to the 
potential energy connected to the size of one particle. Due 
to this rather low kinetic energy, the array of particles is ar- 
ranged in triangular order, except for a few layers at the top. 
Some tests with different values of v lead to similar results as 
long as v is not too large. The larger v, the more particles 
belong to the fluidized part of the system at the top and in 
the fluidized part of the system we can not observe cracks. 
For lower values of v we observed an increasing number of 
events per unit time and thus increasing computational costs 
for the simulation. At t = we remove the bottom, switch 
on dissipation and friction and let the array fall. 

Performing simulations with different initial conditions and 
different sets of material's parameters, we observe strong fluc- 
tuations in position and shape of the cracks. However, the 
intensity or the probability of the cracks seems to depend on 
the material's parameters rather than on the initial condi- 
tions. The behavior of the system depends on friction and on 
dissipation as well: For weak friction we observe only random 
cracks, which would occur even in a dense hard sphere gas 
without any friction, simply due to random fluctuations and 
the internal pressure. With increasing friction, cracks may 
even span the whole system and sometimes be correlated to 
slip planes. Furthermore, we observe cracks to occur more 
frequent for lower dissipation. 

In Fig. ^(a) we present typical snapshots of an experiment 
and of a simulation at time t — 0.04 s. We observe long range 
cracks from both, experiments and simulations as well. For 
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the experiment we use a container of width L/d — 24 with 103 
layers of oxidized aluminum particles. The vertical 2D cell is 
made up of two glass windows for visualisation and of two 
lateral walls made of plexiglass. The gap between the glass 
windows is a little larger than the bead diameter, what leads 
to a small friction between particles and front/back walls, 
while the friction with the side walls may be large. Different 
heights and wall-materials were used and the results could 
be scaled with a characteristic length £ = L/(2Kp w ), with 
the dimensionless parameter K which characterizes the con- 
version of vertical to tangential stresses and the coefficient of 
friction fj, w . The scaling in the regime before cracks occur 
leads to the value Kfj, w « 0.12. For a more detailed descrip- 
tion of the experimental setup and data see Ref. MM. In the 
simulation of Fig. Q we have N — 1562 and the parameters 
L/d = 20.2, e = 0.96, e w = 0.92, fj, = 0.5, fi w = 1.0, and 
A) = Pow =0.2. We varied the coefficients of friction in the 
intervals < fJ, < 1 and < fj, w < 10. Furthermore, we varied 
the coefficients of restitution in the range 0.80 < eore p < 0.98. 
However, the occurrence of cracks is quite independent of the 
parameters used, as long as the coefficients of friction, n and 
/iu, are sufficiently large. Furthermore, cracks occur faster for 
stronger dissipation since a highly dissipative block dissipates 
the initial energy faster. 

For the above parameters, we find - like in experiments - a 
large number of cracks, overlapping and interferring. In the 
upper part of both, experiment and simulation, we sometimes 
observe cracks only on one side. In contrast to Fig. [j](a) we 
present in Fig. [j] (b) a specific simulation with only one strong 
crack, on which we focus in more detail. This crack separates 
the system into a large upper and a small lower part and is 
best visible in Fig. [l](b) at t = 0.068 s. Here we use a reduced 
wall friction, i.e. ^ = fi = 0.5, and stronger dissipation, i.e. 
e w = e = 0.9, while all other parameters, including the initial 
configuration, are the same as for (a). 

In order to distinguish consistently, we will refer to the sim- 
ulation with large simulation (a) , and to the simulation 
with small simulation (b) in the following. 

The important feature of the crack in Fig. |l](b) is that it 
seems to be connected to one single particle. We indicate the 
vertical position of particle #371 with a small bar in Fig. 0(b). 
The crack of simulation (b) is connected to an increase of the 
number of collisions per particle, indicated by the greyscale 
on Fig. |l|(b) . Black or white correspond to zero or more than 
twenty collisions during the last millisecond respectively. Note 
that the increase in the number of collisions is here equiva- 
lent to an increase in pressure. Already at time t = 0.048 s, 
particle #371 peforms more collisions than the average parti- 
cle. The pressure around particle #371 increases and at t = 
0.056 s a region in which the particles perform a large number 
of collisions spans the whole width of the container. We call 
such an array of particles under high pressure 'dynamic arch'. 
At t = 0.060 s the pressure decreases and a large crack opens 
below particle #371. At later times, we observe new pres- 
sure fluctuations in the array. In conclusion, a crack seems 
to begin at one point, i.e. one particle, where the pressure 
increases accidentally. 

Before we look in more detail at the behavior of particle 
#371, we present for convenience, a picture of some selected 
particles around #371 from (b), in Fig. [j](c). 



A. Evidence for Stick-Slip Behavior 

From Fig. [j](b) we evidenced that a crack may originate 
from one particle only. Now we are interested in the velocity 
of one specific particle during its fall. Following the order of 
simulations (a) and (b) we firstly present the case of large wall 
friction (a) and secondly the case of smaller wall friction (b) 
in the following. Remember that particle #371 is the origin 
of one single crack in simulation (b), while its behavior is 
similar to the behavior of many others, close to the boundary, 
in simulation (a). 
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FIG. 2. Plot of the angular velocity, uir, and of the linear, 
vertical velocity, V z , of particle #371 versus time. The sim- 
ulations are the same as in Fig. hl(a) and (b). The line —gt 
corresponds to a freely falling particle. 

Due to gravity, the particle is accelerated downwards and 
after a collision with the right wall it will presumably rotate 
counterclockwise. Therefore, we compare the linear velocity, 
V z , with the rotational velocity of the surface, cor. In Fig. Q we 
plot both, the linear vertical and the angular velocity of par- 
ticle #371 as a function of time. The horizontal line indicates 
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zero velocity and the diagonal line indicates the free fall ve- 
locity, —gt. The full curve gives the (negative) linear velocity 
in vertical direction, V z , while the rotational velocity of the 
surface, lor, is given by the dashed curve. Negative tor values 
correspond to counterclockwise rotation and for V z = lot we 
have the contact point of the particle at rest relative to the 
wall. The particle adapts rotational and linear velocity, or in 
other words, the contact point sticks. This event occurs when 
the two curves merge. In Fig. ^|(a) we observe a small angu- 
lar velocity up to t ~ 0.02 s when particle #371 first sticks. 
Note that in this simulation (a) #371 can not be identified as 
the initiator of one of the numerous cracks; however, it sticks 
and slips several times. An increase in angular velocity goes 
ahead with a decrease of linear velocity due to momentum 
conservation. At larger times we observe the angular velocity 
decreasing, i.e. the particle slips, and short time after, sticks 
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Now we focus on simulation (b) and the behavior of the 
particle, from which the crack started. In Fig. |^(b) we observe 
a small angular velocity up to t ~ 0.050 s. At t — 0.055 s the 
velocities are adapted for some 0.005 s before the particle slips 
again. Since the pressure fluctuations are visible in Fig [j](b) 
already at t = 0.048 s, we conclude that pressure fluctuations 
in the bulk lead to a sticking of a particle surface on the 
wall. This particle is slowed down and thus will perform more 
collisions with those particles coming from above, what leads 
to an increase of pressure. An increase of pressure allows, in 
general, a strong Coulomb friction and thus a sticking of the 
contact point. When pressure decreases, the particle surface 
does not longer stick on the wall. 

Since the sticking might also occur between particles, we 
examine the angular velocity of the particles in the neighbor- 
hood of the sticking particle, i.e. the particles to the upper 
left of particle #371. In Fig. g we plot the angular velocities of 
particles #371, #390, #409, and #428 for the simulations (a) 
and (b). We observe from both figures as a response to fric- 
tion, an auto-organization of the spins as also observed in ID 
experiments and simulations of rotating frictional cylinders 
[ }l2[|l3fl . Spin stands here for the direction of angular veloc- 
ity of a particle. We see that the direct neighbor of #371 to 
the left and upwards, i.e. #390, rotates in the opposite di- 
rection as #371. Thus a counterclockwise rotation of particle 
#371 leads to a clockwise rotation of #390. This is consistent 
with the idea of friction reducing the relative surface veloc- 
ity. The next particle to the upper left, i.e. #409, is again 
rotating couterclockwise, following the same idea of frictional 
coupling. In Fig. ^|(a) we observe only a weak coupling be- 
tween #409 and #428, whereas in Fig. |(b) particle #428 is 
also rotating clockwise, even when the absolute value of the 
angular velocity is smaller. Thus we have not only a stick-slip 
behavior of particles close to the side walls, but also a cou- 
pling of the spins of neighboring particles. If the coupling is 
strong enough, we observe an alternating, but decreasing an- 
gular velocity along a line. We will return to this observation 
in subsection III C . 
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FIG. 3. Plot of the angular velocities of the particles 
#371, #390, #409, and #428 for the same simulations as 
in Fig. bJ(a) and (b). These particles were initially arranged 
on a line, tilted clockwise by 60 degrees from the horizontal, 
sec Fig. 0(c). 



B. Number of Collisions, Kinetic Energy and 
Pressure 

Since the occurrence of a crack is, in general, connected to 
a large number of collisions, we plot in Fig. ^ the number of 
collisions per particle per millisecond, N c , for simulations (a) 
and (b). In Fig. ^(a) we observe an increase in the number of 
collisions, which is related to the first sticking event of particle 
#371. However, particles deep in the array possibly perform 
much more collisions, see #428 in Fig. ^a). Obviously, an 
increase in N c for one particle is connected to an increase 
in N c for the neighbors. In Fig. ^|(b) we find values of N c 
comparable for direct neighbors, i.e. #371 and #390. At t = 
0.056 s we observe a drastic increase in N c which also involves 
the particles deeper in the array, i.e. #409 and #428. 
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FIG. 4. Log-lin plot of the number of collisions, N c , per 
millisecond (ms) as a function of simulation time, t. Plotted 
is N c for the four particles noted in the figure. The data are 
from the same simulations as already presented in Fig. |l| (a) 
and (b). 

In order to illustrate the reduced falling velocity, connected 
to a large number of collisions we plot in Fig. [B| the averaged 
kinetic energy K for particles between heights z and z + dz 
(we use here dz = 0.002 m). We plot K = (l/AQE^w? 
(disregarding the mass of the particles), as a function of the 
height, z. Note that Vi is the velocity of one particle relative to 
the walls, and not relative to the center of mass of the falling 
sandpile. We observe different behavior for the simulations 
(a) and (b): A rather homogeneous deceleration of the pile in 
(a) and one 'dynamic arch' connected to a strong deceleration 
in (b). For strong friction at the walls (a), all particles in 
the system are slowed down due to frequent collisions of the 
particles with the side walls and inside the bulk. For the 
times t = 0.05, 0.06, and 0.08 s we observe, from the bottom 
indicated by the left vertical line in Fig. o, a decreasing K 
with increasing height z up to z m 0.02 m, where the slope 
of K almost vanishes. In Fig. ^(b) the particles are falling 
faster at the beginning until the first dynamic arch occurs at 
t w 0.056 s, see Fig. Mb). The high pressure exerted on the 
walls, together with the tangential friction at the walls, leads 
to a local deceleration, i.e. the dip in the K profile for t — 0.06 
s. This dip identifies a dynamic arch of slow material which 
temporarily blocks the flow. Later in time, the particles from 



above arrive at the slower dynamic arch, which again leads 
to great pressure, a large number of collisions and thus to a 
further reduction of K (see the K profile for t — 0.08 s). 
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FIG. 5. Plot of the kinetic energy K as a, function of 
height z for different times for the same simulations as already 
presented in Fig. |l] (a) and (b). 

In order to understand how the pressure in the bulk is con- 
nected to N c and K we plot in Fig. ^ the particle-particle 
pressure, pp, in arbitrary units as a function of height for the 
simulations (a) and (b). pp is here defined as the sum over 
the absolute normal part of momentum change 



pp{z,t) = Ej AP 



Wl 



(4) 



for each particle in a layer [z, z+dz] in a time interval [t, t—dt]. 
For this plot we use dz — \/3d (what corresponds to a height 
of two particle layers) and the integration time is here dt = 
0.005 s. For strong wall friction and low dissipation (a), we 
observe already at t = 0.02 s a quite strong pressure with 
a maximum close to the bottom of the pile. At t = 0.05 s 
the maximum in pressure moved upwards, not only within 
the falling pile but also in coordinate z. Furthermore, the 
maximum is about one order of magnitude larger than before. 
Later, the pressure peak moved further upwards and at t — 
0.06 s begins to decrease until at larger times the array is 
dilute and pressure almost vanished. For weak wall friction 
(b) we find only a weak pressure at time t = 0.02 s until at 
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time t = 0.06 s a sudden increase of pressure, connected to the 
increase in N c and the decrease in K, appears. The occurrence 
of two pressure peaks of different amplitude (small pressure 
at the bottom and large pressure at the top) is consistent 
with the predictions of Duran et al. p^ , which state that a 
small lower pile will be decelerated less than a large upper 
pile. Thus the crack continues opening. 
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pendicular to a line of strong pressure, such that the order in 
Fig. (?](a) indicates an arch like structure. 

In Fig. 0(b) we plot snapshots from simulation (b) and 
plot the change of velocity, i.e. Av x — v x (t + St) — v x (t) and 
Av z = v z (t + St) — v z (t) + g St with St — 10 -3 s. Comparing 
Figs. |l|(b) and Q(b) we identify the large number of collisions, 
starting from particle #371, with a momentum wave propa- 
gating from #371 towards the left wall and also diagonally 
upwards. When the momentum wave arrives at the left wall 
it is reflected and moves mainly upwards. We relate this to 
the fact that the material below the dynamic arch is falling 
faster than the dynamic arch, such that not much momentum 
change takes place downwards. After several milliseconds the 
momentum wave is not longer limited to some particles only, 
but has spread and builds now an active region with great 
pressure, i.e. the dynamic arch. 
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FIG. 6. Log-lin plot of the pressure, pp, as a function of 
height z for different times for the same simulations as already 
presented in Fig. m (a) and (b). 



C. Long Range Rotational Order and Momentum 
Waves 

We have learned from Figs. ^ and |^ that, connected to a 
large number of collisions, the direction of the angular veloc- 
ity, i.e. the spin of the particles, may be locally arranged in 
an alternating order along lines of large pressure. In order to 
proof that this is not only a random event we plot in Fig. [j](a) 
some snapshots from simulation (a) and indicate clockwise 
and counterclockwise rotation with black and white circles 
respectively. We observe, at least in some parts of the sys- 
tem, that spins of the same direction are arranged along lines. 
The spins of two neighboring lines have different directions. 
The elongation of the ordered regions may be comparable to 
the size of the system. Note, that lines of equal spin are per- 



(a) t = 0.04 s 
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t = 0.053 s 0.054 s 0.055 s 0.056 s 0.057 s (b) 



FIG. 7. (a) Snapshots of the simulation from Fig. |l|(a) at 
different times. The greyscale indicates here the direction of 
the angular velocity, i.e. black and grey correspond to clock- 
wise and counterclockwise rotation respectively, (b) Snap- 
shots from the simulation from Fig. |l|(b) at different times. 
The lines indicate for each particle the change in velocity due 
to collisions within the last millisecond, i.e. the last image. 

Examining the rotational order in simulation (b), we ob- 
serve that the spin order is a consequence of the momentum 
wave. In our model, strong coupling is related to a large num- 
ber of collisions and thus to a large pressure. This is due to 
the fact, that informations about the state of the particles are 
exchanged only on contact. Therefore, lines of equal spin are 
mostly perpendicular to the lines of great pressure, a finding 
that is also discussed in Ref. Il4|. 



IV. SUMMARY AND CONCLUSION 

We presented simulations of a 2D granular model mate- 
rial, falling inside a vertical container with parallel walls. We 
observe fractures in the material which were described as a 



so-called 'discontinuous decompaction' which is the result of 
many cracks breaking the granular assembly into pieces from 
the bottom to the top. The case of high friction and quite 
low dissipation (a) is a system which shows the same behav- 
ior as also found in various experiments. Many cracks occur 
and interfere. When we investigate an exemplary situation 
with rather low friction and high dissipation (b) we observe 
one isolated crack. We followed in detail the events which 
lead to this crack. Due to fluctuations of pressure (equiv- 
alent to fluctuations in the number of collisions, N c ), some 
particles may transfer a part of their linear momentum into 
rotational momentum. This happens when the surface of a 
particle with small angular velocity sticks on the wall. Stick- 
ing means here that the velocites of particle surface and wall 
surface are adapted. 

The momentum wave, starting from such a particle, leads 
to a region of large pressure, which spans the width of the 
system, i.e. a dynamic arch. Due to the strong pressure, the 
dynamic arch is slowed down by friction with the walls. The 
material coming from above hits the dynamic arch such that 
a density and pressure wave, propagates upwards inside the 
system. 

Under conditions with quite strong wall friction and rather 
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weak dissipation, the fluctuations in the system and also the 
coupling with the walls are greater, such that many particles 
are sticking on the walls. This leads to several pressure waves, 
interferring with each other such that the system is slowed 
down more homogeneously. 

With the dynamic simulations used, we were able to re- 
produce the experimentally observed discontinuous decom- 
paction and to propose an explanation, how the cracks 
occur. The open question remains, if the situation discussed 
here is relevant for all types of experiments in restricted ge- 
ometries. The features described here, i.e. rotational or- 
der, stick-slip behavior, and momentum waves are not yet 
observed in experiments. Besides, the phenomenon of stress 
fluctuations has been shown to be important for both, the be- 
havior of static |H| and of quasistatic granular systems JlBj ]. 
Furthermore, the behavior of cracks in polydisperse or three 
dimensional systems is still an open problem. 
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